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I. INTRODUCTION 



The iV(1535) and iV(1650) resonances appear as outstanding features not only in elastic 
ttN scattering in the strangeness zero Su (L 2 t,2j) partial wave but also in other meson- 
baryon reactions at intermediate energies. In quark model approaches these excited nucleon 
resonances are mainly composites of three valence quarks, and their widths are computed 
as matrix elements of hadronic transition operators. However, the description of hadron 
scattering reactions becomes cumbersome in this framework, and it requires quite elaborated 
techniques as the resonating group approach, where it becomes extremely difficult to impose 
Chiral Symmetry (CS) [§. 

Renouncing to find out a picture of the hadron as a valence quark bound state, a different 
point of view consists of describing scattering reactions taking the hadrons as the relevant 
degrees of freedom at low energies. Then, resonances manifest themselves as poles of the 
scattering amplitude in a certain Riemann sheet in the complex energy plane. To perform 
such a program requires to implement unitarity in the model. A multichannel K— matrix 
method is used in the work of Ref . |J . Though CS is not incorporated, this phenomenological 
approach is able to reproduce a large amount of data related to the nN — > ttN reaction. 
The, three body final state, two pion production channel ttttN is incorporated through an 
effective use of two body channels with higher mesonic and baryonic resonances. In this 
paper, we will work in this latter type of approaches, but explicitly imposing CS constraints 
as an indirect way of incorporating the bulk of the underlying Quantum Chromo Dynamics 
(QCD). Thus, we will establish a unitarity scheme based on the Chiral Perturbation Theory 
(ChPT) amplitudes. 

CS provides important constraints to the description of low energy hadronic processes 
and, in particular, to baryon-meson dynamics. There have been previous studies of the 
ttN— Su partial wave using a coupled channel formalism and imposing CS constraints. In 
Ref. || a Schrodinger coupled channel treatment was employed and the additional inclusion 
of phenomenological hadronic form factors was invoked. Within this framework, the p— wave 
contribution has also been recently examined |4|]. The two pion production channel is not 
considered in these works. Perturbative estimates [|5|- |T] for the reaction ttN — > TnvN 
indicate that this three-body channel keeps moderately small not only at threshold but also 
in a certain region above it. An attempt to include the two pion production reaction can be 
found in Ref. j8|, but the treatment of the tttcN channel is effective and it is represented by 
a unphysical two body channel which represents all remaining processes. 

In Ref. U the Bethe-Salpeter Equation (BSE) has been employed in the spirit of an 
Effective Field Theory (EFT). There, the rnrN channel is not considered either and the 
authors require the introduction of less renormalization constants than allowed by CS. De- 
spite of these restrictions, the model describes not only the elastic ttN channel, but also the 
two-body inelastic ones in an energy window around the iV(1535) resonance. Nevertheless, 



it fails at threshold ||10|| . Given this partial success and the great technical difficulties to 
solve the BSE incorporating the three body ttkN channel, one might wonder what features 
of the data can be explained incorporating CS constraints and restoring two body unitarity. 

In the present work, we restrict our study to the non-strange meson-baryon Su partial 
wave and adopt a similar framework as in those references, but with some important differ- 
ences. First, we will implement exact unitarity by solving the BSE taking the needed input 
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from lowest order relativistic ChPT. A similar program has been successfully undertaken 
both in the pion-pion sector ['J 1] and in the ttN—P 33 partial wave [0]. Thus, we avoid the 
use of phenomenological form-factors and all required information (low energy constants) 
can be, in principle, obtained from higher orders in the chiral expansion. Besides, we aim 
at describing not only a narrow energy window, placed at threshold or in the neighborhood 
of some resonance, but also a wider energy region ranging from tcN threshold up to almost 
a Center of Mass (CM) meson-baryon energy of ^/s = 2 GeV. 

As we discussed at length in Ref. [nj the BSE, in the context of EFT's, can be solved 



in two different schemes: off-shell and on-shell. Here, we use the off-shell scheme because 
of the lack of information on the next-to-leading order in the chiral expansion. In this 
scheme, the on-shell scattering amplitude requires some knowledge of the off-shell behavior 
of the two particle irreducible amplitude (potential) . After renormalization of the amplitude 
this off-shell input leads to a finite number of phenomenological constants which encode the 
detailed underlying short-distance dynamics. In practice, these constants can be either fitted 
to experiment or determined by matching the resulting Bethe-Salpeter (BS) amplitude to 
standard ChPTQ. Obviously, the method of determining the constants by matching to ChPT 
seems a better one than a direct fit to experimental dataQ. For the case of meson-baryon 
scattering, the only known information coming from ChPT involves tree-level amplitudes 
and free propagators; there is no possibility to compare with ChPT beyond leading order 
and thus one is forced to fit the unknown low energy constants (LEC's) to data. 

As it is the case in the purely mesonic sector, the off-shell scheme generates a rich struc- 
ture of unknown constants which allow for a good description of the amplitudes. Although 
the appearance of more undetermined constants may appear a less predictive approach as, 
say, putting a cut-off (one single parameter) in the divergent integrals as it is done in Ref. , 
it reflects the real state of the art of our lack of knowledge on underlying QCD dynamics. 
The number of adjustable LEC's should not be smaller than those allowed by symmetry; 
this is the only way both to falsify all possible theories embodying the same symmetry prin- 
ciples and to make wider the energy window which is being described. Limiting such a rich 
structure allowed by CS results in a poor description of experimental data. 

Before going further we would like also to say a word on the opposite situation, i.e., 
the possibility of having more low energy parameters than one actually needs. A possible 
redundancy of parameters is obviously a undesirable situation. In standard ChPT the 
number of LEC's is controlled to any order of the calculation by crossing and unitarity. 
Moreover, the dependence of the observables on them is strictly linear, so that it becomes 



1 Ideally, these phenomenological parameters should be computed from first principles, a yet im- 
possible task. 

2 In addition, if the matching is possible the off-shell scheme becomes unnecessarily complicated, 
as compared to other methods directly unitarizing the final on-shell amplitude given in terms of 
the standard low energy constants. For details, on the on-shell BSE approach or on the Inverse 
Amplitude Method (IAM) for meson-meson scattering, see for instance, the thorough discussion in 
the second entry of Ref. Jll|. For ttN elastic scattering, recent IAM studies have been pursued in 
Refs. @ and 0. 
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possible to detect such a redundant combination in case it occurs^. In a unitarized approach, 
the only way to avoid this parameter redundancy is to match the unitarized amplitude to 
the standard ChPT amplitude. As we have already said, there is no standard one loop 
ChPT calculation for the Su partial wave of meson-baryon reaction with open channels to 
compare with. An indirect way to detect such a parameter redundancy might be through a 
fit to experimental data if the errors in some parameters turn out to be very large. 

We have considered four coupled channels: ttN, r]N, KA, KT, and taken into account 
SU(3)— breaking symmetry effects but neglected the considerably smaller isospin violation 
ones. We have found that CS allows for a solution of the BSE which is flexible enough 
to describe the elastic ttN phase-shifts, and the ix~p — > rjn and — ► K°A cross section 
data from threshold to CM energies well above the iV(1650) resonance. Besides the rest of 
elastic and inelastic two-body reaction channels^ implicit in the adopted formalism come 
as predictions of the model. The position of the complex poles in the second Riemann 
sheet of the amplitudes determine masses and widths of the Sii— iV(1535) and — iV(1650) 
resonances, which turn out to be in reasonable agreement with experiment. Preliminary 
results were presented in | fL8|1 . 

The paper is organized as follows: In Sect. we present the basic formalism used along 
the paper. We start with the chiral Lagrangian relevant to our calculation, from which 
the lowest order meson-baryon two particle irreducible matrix amplitude is deduced. After 
presenting our notations for the coupled channel kinematics we discuss the basic pertinent 
features of the BSE for s— wave meson-baryon scattering. Using the amplitude from lowest 
order ChPT as the potential, we solve and renormalize the BSE in the spirit of an EFT. 
In Sect. [Ill] we present our numerical results, together with a detailed discussion on the 
fitting procedure and Monte-Carlo propagation of inherited error bars to all possible reaction 
channels. The quantum field theoretical interpretation of resonances as unstable particles 
requires determining their mass and width as poles in a unphysical Riemann sheet. In 
our case there are 16 sheets which we discuss in some detail, and we search for the most 
important pole singularities. Error estimates are made in terms of the available experimental 
uncertainties in the phase-shifts and amplitudes. Finally, in Sect. IV we present some 
conclusions and outlook for future work. 



3 A good example of this is tttt scattering to two loops ( Jl5[ , (Tq]) where one gets, besides the 
four one loop parameters £1,2,3,4, six new parameters but in redundant combinations. Instead, it 
is customary to use the six 61,2,3,4,5,6 independent combinations, which depend on the one loop I's 
and the six new two loop parameters and contain mixed orders. Such a situation also takes place 



in irN scattering in HBChPT at fourth order [17 



4 Each of the entries of the 4x4 matrix solution of the BSE is the T— scattering amplitude for a 
meson-baryon reaction constructed out of the four considered channels: ttN — * irN, ttN — * r]N , 
irN -> KA, irN -> KT,, t)N -> 77JV, r)N -> KA, t]N -» KT,, KA -► KA, KA -► KE, KE -» KE 
and the reverse processes. 
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II. THEORETICAL FRAMEWORK 



A. Chiral Baryon-Meson Lagrangian 

At lowest order in the chiral expansion the chiral baryon meson Lagrangian contains 
kinetic and mass baryon pieces and meson-baryon interaction terms and is given by 



C x = Tr [B (if - M B ) b) + \ V Tr {fl 7 " 75 {u„ B}} + ^Tr {B^ l5 [u„ B]} , (1) 



The meson kinetic and mass pieces and the baryon mass chiral corrections are second order 
and read 



£ 2 = ^Tr{^ + (C/t x + x t[/)} 

- b Tr( X +)Tr(BB) - ^Tr (B X +B) - b 2 Ti(BB X , 
where "Tr" stands for the trace in SU(3). In addition, 

V^B = d^B + 1 [ u%u + ud^ , B } , 



(2) 



U = u 2 = e iV ~ 2 */ f 



X+ = u ] xu ] + ux ] u . 



X = 2B M 



(3) 



Mb is the common mass of the baryon octect, due to spontaneous chiral symmetry breaking 
for massless quarks. The SU(3) coupling constants which are determined by semileptonic 
decays of hyperons are T ~ 0.46, V ~ 0.79 [T + V = g& — 1.25). The constants B 
and / (pion weak decay constant in the chiral limit) are not determined by the symmetry. 
The current quark mass matrix is M. = Diag(m u , m^, m s ). The parameters bo, b\ and 62 
are coupling constants with dimension of an inverse mass. The values of &i and 62 can 
be determined from baryon mass splittings, whereas 60 gives an overall contribution to the 
octect baryon mass Mb- The SU{3) matrices for the meson and the baryon octect are 
written in terms of the meson and baryon spinor fields respectively and are given byQ 



7T~ 

K~ 



75^ 



7T 

K c 



and 



B 



. J_V0 _j_ J_ A 

wo 



K+ \ 
K° 

-fa/ 



P 
n 

7E A 



(4) 



(5) 



5 For the purpose of our work we do not consider any mixing between octect and singlet SU(3) 
representations 
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respectively. The MB — > MB vertex obtained from the former Lagrangian readsQ 

Cmb^mb = ^Tr [BY [ [ $ , ] , B ] } . (6) 

Assuming isospin conservation, the scattering amplitude-] in the Dirac spinor basis, which 
relation to the cross section is given in the next subsection, at lowest order is given by 

t?(k,k') = ^(k' + k") (7) 

where k and k' are incoming and outgoing meson momenta and D a coupled-channel matrix. 
For strangeness S = and isospin T = 1/2 the coupled channel matrix D reads[| 



D 



T=l/2 
5=0 



irN 
-2 


-3/2 
V+l/2 



r)N 




+3/2 
+3/2 



KA 
-3/2 
-3/2 







KE 

+ l/2\ 
+3/2 


-2 / 



TjN 

KA 
KY, 



(8) 



While amplitudes follow the chiral symmetry breaking pattern from the effective La- 
grangian to a good approximation, it is well known that physical mass splittings have an 
important influence when calculating the reaction phase space. Indeed, the correct loca- 
tion of reaction thresholds requires taking physical masses for the corresponding reaction 
channels. We have taken into account this effect in our numerical calculation. Besides, 
chiral corrections to the amplitudes also incorporate explicit CS breaking effects in addition 
to those already present in the lagrangians above. In practice, we use different numerical 
values for f n , fx and f v , as it is discussed in Sect. |T|. This can be easily accomplished 
through the prescription 



D/f - r i D}- 



f = Diag (f n , f v , f K , f K ) 



(9) 



For simplicity and a more clear book-keeping of chiral order dependences we will use the 
D/f 2 notation throughout the paper, meaning Eq. (|9|) in practice. 



We have omitted the pieces proportional to the couplings V and T because they do not contribute 
to s— wave. On the other hand, the lagrangian below does not lead to a pure s— wave contribution 
and a further projection will be required. 

7 We use the convention, in symbolic notation, —iTmb-^mb = +i£MB-+MB- 

8 There is a mistake in the relative phases of Ref. Q. We thank A. Ramos for confirming this 
point to us. We use the isospin phase convention of Ref. [pcf] : negative phases for the isospin states 
-|tt+> , -\K ) , -|E+), -IH 3 ) , -\TF) , -\n). 
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B. Scattering Amplitude and Kinematics 



The coupled channel scattering amplitude for the baryon-meson process in the isospin 



channel, T — | 



B(M A , P-k, s A ) + M(m A , k) -> B(M B , P - k', s B ) + M(m B , k') 



(10) 



with baryon (meson) masses M A and M B ( m A and m B ) and spin indices (helicity, covariant 
spin, etc..) s A , s B , is given by 



T P [B{k', s B } <- A{k, s A }] = u B (P - k', s B )t P (k, k')u A (P - k, s A ) 



(11) 



Here, u A (P — k, s A ) and u B (P— k', s B ) are baryon Dirac spinors[] for the ingoing and outgoing 
baryons respectively, P is the conserved total four momentum and tp(k, k') is a matrix in the 
Dirac and coupled channel spaces. On the mass shell and using the equations of motion for 
the free Dirac spinors (f—fk — M A )u A (P — k) = and its transposed u A (P — k) (f — k 1 — M A ) = 
the parity and Lorentz invariant amplitude t P can be written as: 



t P (k, A;') |on-sheii = h(s,t)f + h(s,t) 



(12) 



with s = P 2 = f 2 , t = (k — k') 2 and t\ and £2 matrices in the coupled channel space. The 
normalization of the amplitude Tp is determined by its relation to the CM differential cross 
section, and it is given above threshold, s > Ma.x{(M A + m A ) 2 , (M B + m B ) 2 }, by 



^ [B{k B , s B } <- A{k A , s A }} = 1 

dil QAtt 2 s \k A \ 



T P [B{k B , s B } <- A{k A , s A }} 



(13) 



Rotational, parity and time reversal invariances ensure for the on shell particles 



T P [{k B , s B } <- {k A , s A }] = -8n^\ l p\ U(s, 6)5 SASB + i B(s, 9) (n ■ a) SBSA } (14) 

A and B are matrices in the coupled channel space, 9 the CM angle between the initial 
and final meson three momenta and n a unit three- vector orthogonal to k A and k B . Partial 
waves (matrices in the coupled channel space), f[(s), are related to A, B by |21| 



a( S ,9)=y: 

L 



Pl(cos9) 



f. 



fl 



dPh (cos 1 
d9~ 



(15) 



In terms of the matrices t\ and £2 defined in Eq. (|i2[) , the s— wave coupled-channel matrix, 
1 

/o 2 ( s ), is S iven bv: 



9 We use the normalization uu = 2M. 
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/«?(«) 



57TV S 



\k 



B\ 



\k f 



E b + M b JEa + M a 



l r l 



dcOS0 (y/sti(s : t) + t2(-S,t) 



BA 



(16) 



where the CM three-momentum moduli read 



2^~s 



i = A,B 



(17) 



with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz and Ea,b the baryon CM energies. The 
phase of the matrix Tp is such that the relation between the diagonal elements (A = B) in 

the coupled channel space of / 2 (s) and the inelasticities (77) and phase-shifts (5) is the usual 
one, 









- ,-,( 




aa 2i\k A \ 



VA{s)e 



2i6 A (s) 



(18) 



Hence, the optical theorem reads, for s > (Ma + it-ia) 



An 



Im 



/<?(*) 



= °B^A = ^Y. 
AA B B 



/«?(*) 



BA 



<J AA + 



71 



(19) 



where in the right hand side only open channels contribute. 



C. Bethe-Salpeter Equation 



The Bethe-Salpeter equation reads 



t P (k,k') =v P {k,k') +i 



d 4 q 



tp(q,k')A(q)S(P-q)v P (k,q) 



(20) 



where tp(k, k') is the scattering amplitude defined in Eq. (|TT|) , and vp(k, k') the two particle 
irreducible Green's function (or potential ), and S(P — q) and A(g) the baryon and meson 
exact propagators respectively. The above equation turns out to be a matrix one, both in 
the coupled channel and Dirac spaces. The resulting scattering amplitude tp(k,k') fulfills 
the coupled channel unitarity condition 



xt P (q,k') 5 H 



t P (k,k') -t P (k',k) = 



2 ~ 2 

q — m 



-i(2rr) 
(P-qf 



d 4 q 
- M 2 



t P (q, k) 



(21) 



where ip(k,p) = 7o^p(&,p)7o an d t P (k,p) stands for the total adjoint in the Dirac and 
coupled channel spaces (including also the change s + ie — > s — ie) and rh and M the meson 
and baryon (diagonal) mass matrices respectively. Finally, 5 + (p 2 — m 2 ) = Q(p°)5(p 2 — m 2 ), 
being 9 the Heaviside step function. 
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If the on shell amplitude depends only on the total momentum^ P, tp(k, fc')lon-sheii = 
t(f), as it will be the case below, the unitarity condition can be rewritten in a much simpler 
and useful form as a discontinuity equation above the corresponding physical thresholds 

Disc[t(f = -Disc[J(f )] with DiscL4(s)] = A(s + ie) - A(s - ie) (22) 

where the quadratic and logarithmically divergent integral 

Jir) = i I (trf-^p-l-M (23) 

has been introduced and Cutkosky's rules used to evaluate its discontinuity This integral 
is treated in detail in Appendix |A]. As usual, we take the ie prescription m 2 — > m 2 — ie and 
M — > M — ie which we implicitly assume in the sequel. 

D. Solution of the BSE Equation at lowest order 

The BSE requires some input potential and baryon and meson propagators to be solved. 
We proceed here along the lines proposed in previous work |TTJ and use a chiral expansion 



to determine both potential and propagator. From the chiral Lagrangian one gets at lowest 
order (Eq. fl7|)) 

Vp (k,k')=t?(k,k') = j 2 (k' + k") (24) 

with D the coupled-channel matrix, which was given in Eq. (H). The propagators at lowest 
order are simply the free ones, 

A(q) = - 1 - , S(P -q) = 7 - (25) 

W q 2 -m 2 K H! f-fr-M K ' 

which are diagonal in the coupled channel space. Once we have approximate expressions 
for both the potential and the meson and baryon propagators, we proceed to exactly solve 
the BSE. Second order Born approximation to the solution of the BSE (i.e. approximating 
tp(q, k') in the kernel of the equation by vp(q, k') ) suggests the following form for the exact 
solution: 

t P (k, k') = a{P ) + ffb R (f) + b L (P)^ + tfc^k 1 (26) 

where a,bp,bL and c are Lorentz scalar matrices in the Dirac and the coupled channel spaces. 
At this lowest order of the BSE approach, these matrices only depend on f, thus they turn 
out to be independent of the Mandelstam variable t. 
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It is to say, the functions ti and £2 in Eq- fll^) do not depend on the Mandelstam variable t. 



S 



k k' 

k k ' 

\ q / 

\ / + \/~\/ + \/' V V + . . . 

it u a I E 

P-k p_k' P-k P-k' 



FIG. 1. Diagrams summed by the Bethe Salpeter equation at lowest order. Kinematics denned in the 
main text. 

Such an amplitude contains an infinite sum of diagrams, as shown in Fig. p], but does 
not contain all possible one loop dependences, for instance those coming from vertex renor- 
malization. On the mass shell we may set $ — > f — M and \f — > f — M. Hence, thanks to 
the equations of motion the on-shell amplitude becomes a function of the total momentum 
f and reads 

t(tf) = a + (f - M)b R + b L (f - M) + (f - M) c(f - M) (27) 

where the explicit dependence on If of the a, b^, c matrix functions has been suppressed 
for simplicity. 

Plugging the ansatz for the off-shell amplitude, t P (k,k'), as given in Eq. ( PS| ) into the 
BSE equation and after some algebraic manipulations described in detail in the Appendix ||, 
we get for the inverse on-shell amplitude 

t{f r 1 = - j(f ) + + a- 1 

A = j 2 {f- M, D} + + l(f - M)D-^D(r - M) = »(s) + v(s)f 

= 4{ - [D, M}^^-[D, M] - f 2 {M, D} + - DA^DM + [D, M]A a d} 
j s — M z ' 

ji{-[M, D}-^- 2 [M, D] + DAmD + 2fD) (28) 



vis) 



As we can see, the on-shell unitarity condition expressed in Eq.(^) is manifestly fulfilled. 
To proceed further we can decompose the inverse on-shell amplitude in the form 

t{f)- 1 ^K 1 {s)f + K 2 (s) (29) 

where Ki(s) and K 2 (s) are matrices in the coupled channel space. Straightforward calcula- 
tion yields, 

T r / \ s - rh 2 + M 2 _ . . A A A . . -, . -, 

K l{ s) = J Q (s) + — ^ - -*± + (us - r V) 

K 2 (s) = -MJq(s) + MA " - {vs\T x v - fi)- 1 (30) 

s — M 2 
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Thus, the amplitude can be written in the form of Eq. fll2|) 



t(f)=t 1 (s)f + t 2 (s) 



(31) 



with 



h = (K lS - K 2 K~ l K 2 ) 
t 2 = (K 2 - K x K^K x s) 



(32) 



As we already said, when the two particle irreducible amplitude,^, and the meson and 
baryon propagators entering in the BSE, the lowest ChPT order ones are taken, the matrices 
(in the coupled channel spaces) t\ and t 2 turn out to be independent of the t— Mandelstam 
variable^. Hence the angle integral in Eq. (|l^) becomes trivial and apart from kinematical 
factors, the relevant combination entering in the s— wave scattering amplitude is 



t(s) = ti(s)Vs + t 2 (s) = (^i(s)Vs + K 2 {s] 



(33) 



After some algebraic reshuffling the expression for the inverse coupled channel matrix am- 
plitude can conveniently be written as 



Us 



(VS + M) : 



m 



-■Us) + 



A 



riiAI 



+ 



2^ ' y/s-M 



2^ 



-D 



i 

Finally, the s— wave coupled-channel matrix amplitude / 2 (s) reads 



M 



(34) 



fHs) 



1 



BA 



jTT V s 



E b + MbJE a + M a 



t(s) 



BA 



(35) 



At this point we have to renormalize the divergent integrals appearing in the solution of 
the BSE. This issue has been carefully discussed in the context of meson-meson scattering 
(see Sect. 3.4 of Ref. |ll[) and applies equally well in the present context. In summary, 
the result of that discussion amounts to write the renormalized amplitudes (finite) in the 
same way as the divergent amplitudes but with the previously divergent integrals taken as 
finite renormalized constants. Ideally one would like to derive them from the underlying 
QCD dynamics, but in practice it proves easier to fit them to experiment. This amounts to 
consider, besides the physical masses and weak meson decay constants 12 fitting parameters 
that define three diagonal matrices in the coupled channel space which appeared already in 
the solution of the BSE given, e.g., in Eq. (E3). These matrices are 



11 Despite of that and because the Dirac structure the amplitude of Eq. ([H]) not only contains 
s— wave, but also a small p— wave. 
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J (s = (m + Mf 



/ 7 


n n 
u u 




n 
U 


t n 

J V N U 


n 
U 


U 


U Jka 


U 


V 





u K 2j / 


/ A 

( ^N,l 


n n 
u u 


n \ 


U 


A n 


U 


U 


n A 
U A A 


U 


i n 


n n 

u u 


/ 


(K 





\ 





A„ 








A K ,i 





V o 








A T 



where we have denoted the meson-baryon low energy constants Jq(s 
ir,r],K,K and j = N, N, A, S of Eq. (K3) as J«. The matrix A 



(36) 



(mi + Mj) 2 ), % = 
which appears in 
^. On 



Eq. (p9|) is determined by the matrices A„ and A^ above, as it is defined in Eq. ([A 
the other hand, the A^i and the Ak,2 matrix elements of the matrix A,^ and the Ajv,i and 
the Ajv,2 matrix elements of the matrix A^ have been taken in general to be different. This 



is because, though their formal expression as divergent integrals in Eqs. (A5) and (|A7|) are 
the same, after doing the needed renormalization there is no reason why the finite parts in 
these two channels should coincide^- 



III. NUMERICAL RESULTS 

Throughout the paper we will use the following numerical values for masses and weak 
decay constants of pseudoescalar mesons (all in MeV), 

= 547.45 m K = 497.67 
M A = 1115.68 M s = 1192.55 
U = f K = 1.3/. (37) 

where for the channel 11 the proton mass is used, because the data have been obtained from 
the 7i~p scattering, and for the channel 22 we take the neutron mass, because the available 
data come from ir~p — > rjn. In this way we ensure the exact and physical position of the 
thresholds. This proves important due to the proximity of the A^(1535) to the rjn threshold. 



M, 



N 



rriir = 139.57 
M p = 938.27 



Mr, 



= 939.57 
A = 93.2 



12 This point is clearly exemplified in the irir BSE treatment, see Eq. (A. 15) in Ref. [11], where 
constants which stem from the same divergent integrals, after renormalization become in fact 
different functions of the SU(2) low energy constants Vs. 
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A. Fitting procedure 

We perform a x 2— fit with 12 free parameters considering the following experimental data 
and conditions: 



S u 71 N elastic phase shifts and inelasticities p2|, 1077.84 MeV < yg < 1946.52 MeV: 



In this CM energy region, there are a total number of 281 phase shifts and inelasticity 
data points . Though we have considered four coupled channels, the three-body rnvN 
channel is not explicitly considered. This omission influences both the phase shifts 
and the inelasticities and we will assume here that the effect is much more important 
for the inelasticities than for the phase shifts. Thus, we have fitted the phase-shifts 
while inelasticities have been considered only to impose some constraints on the fit. 
While 7] > 0.99, we have considered that the ttttN channel is essentially closed. In the 
data, this is the case for CM energies below a/s = 1406.4 MeV. In this energy region, 
we have assigned to the phase shifts a 3% relative error added in quadrature with a 
systematic 1° absolute error, in the spirit of Ref. ||. In this way, we are assuming that 
any ttttN subthreshold effects are effectively incorporated in the systematic error men- 
tioned above. At higher energies, inelasticities are smaller than 0.99 and we provide 
the phase shifts with a systematic 15° absolute error added in quadrature with a 3% 
relative error. The reason for this big systematic error is to account for the explicit 
omission of the, now open and likely important, three-body channel. 

Despite of being able to account for a part of the inelasticities (rjN, KA and A"S 
channels ) and because of the explicit omission of the ttttN channel, inelasticities have 
not been fitted. Nevertheless, some constraints are imposed in order to prevent the 
occurrence of smaller inelasticities than the experimental ones. On a quantitative level 
this means the following. Firstly, we provide the inelasticities with a 3% relative error 
added in quadrature with 0.01 absolute error. In the second step within the x 2 — fit 
procedure and given a set of parameters, we compute the theoretical inelasticity for 
each a/s. If it turns out that for one CM energy the inelasticity is smaller than the 
experimental value, taking into account the provided errors, we strongly disfavor this 
set of parameters by decreasing the total error on the inelasticity for this CM energy 
by an order of magnitude when calculating its contribution to the total x 2 . Besides, 
those energies for which the theoretical inelasticities fall above the experimental ones 
are set to contribute zero to the total \ 2 - I n this way, we do not force the fit to pass 
through the experimental inelasticities at all, but avoid the unphysical scenario where 

^theoretical - experimental 
2^i= v N,KA,KZ °i > "inel 



Total n-p -> rjn cross section g3|, 1488.4 MeV < y/s < 1563.8 MeV 



We fit the region close to the rjn threshold and in terms of the commonly used (/lab 
(incoming pion momentum in the laboratory (LAB) system), the above range corre- 
sponds to 687 MeV ^ ^lab ~ 812 MeV. There is a total number of 11 data points. 



The experimental uncertainties are provided in Ref. ||23|| . In addition, the experimental 
cross section has the contribution not only of the s— wave, object of this work, but 
also of the rest of higher partial waves. Next to threshold the s— wave is the domi- 
nant contribution and the higher energy cut (1563.8 MeV) determines the region up 
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to where it is still a good approximation to the total cross section. For higher energies 
the p— wave does play an important role and cannot be neglected, Ref. [Q. 

We have neglected any possible effect stemming from the tttcN intermediate state in 
this inelastic channel, as it has been also assumed previously in Refs. [§] and |J. 

Total Tr-p -» K°A cross section gfl, 1617.5 MeV < < 1724.8 MeV . 

We fit the region close to the K°A threshold, 904 MeV < g LAB % 1097 MeV with the 
experimental error bars provided in Ref. 3^] . There is a total number of 45 data points 



and the remarks concerning both the contribution of the tcttN channel and p— wave 
effects of the previous item apply also here. 

Note that we have not fitted the ix~p — > K°TP total cross section because of the likely 
sizable isospin 3/2 contribution. 
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FIG. 2. Su— elastic irN phase-shifts and inelasticities as a function of CM energy y/s. Data from 
Ref. p^ ]. Solid lines stand for the lowest-order BSE results with parameters given in Appendix |c[ Dot- 
ted-dashed vertical lines in the bottom plot indicate the energies for which new channels are opened. 
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FIG. 3. Total ir~p — ► r/n and ir~p — > K°A cross sections as a function of the CM energy. Data from 
Rcf. 1 23 1 . Solid lines stand for the lowest-order BSE Su— results with parameters given in Appendix |(]. 
Data for energies above the vertical dotted-dashed lines have not been included in the fit. 



B. Results of the best \ 2 ~ fit 

The best fit parameters are compiled in Appendix The errors on the parameters 
turn out to be fairly small. They are purely statistical and have been obtained from the 
68% confidence level on the best-fit parameter 12- dimensional distributions. We generate 
these parameter distributions out of iV = 10 4 samples. Each of the samples is obtained 
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from a x 2 — fit to a synthetic set of data-points which are obtained from the original one 
by a Gaussian sampling procedure, i.e. the total number of fits is iV = 10 4 . We use the 
so obtained distributions to evaluate the correlation matrix between the parameters. The 
correlation matrix is also given in Appendix 0. Besides, quantities coming from a % 2 — fit, 
are Gauss distributed in the limit of small errors, as it seems to be the case here, therefore 
the correlation matrix determines the parameter distributions. 

Systematic errors on the parameters are not included and they are difficult to estimate. 
We can not completely discard that they might be sizeable. This situation reflects the 
present status of the art in unitarized calculations. 

In Figs. and we show the results of our approach, with parameters given in Ap- 
pendix 0, for those quantities which have been fitted to. In Fig. ||] only the data for energies 
below the vertical dotted-dashed lines have been included in the fit. The overall descrip- 
tion is remarkably good and that gives us some confidence on the used non-perturbative 
resummation procedure based on the BSE. For the elastic ttN — > ttN scattering length we 
get 



fo 2 (s = (m n + M N ) 2 ) 



= 0.179 ±0.004 fm (38) 

ttN^ttN 



where the error is statistical and it has been obtained from those in the best fit parameters 
(Eq. ( |Cip ), taking into account the existing statistical correlations, through a Monte-Carlo 
simulation. This value should be compared both to the recent experimental one 0.252±0.006 



fm of Ref. pj and to the HBChPT result to third order of Ref. |TJ 0.19 ± 0.05 fm. The 
agreement between our coupled channel unitarized scattering length with that from NNLO- 
HBChPT is satisfactory from a theoretical viewpoint since in both cases the same set of low 
energy data have been used. The discrepancy of our number with the experimental one 
of Ref. [24| possibly points toward a too conservative error assignment of the low energy 
phase-shifts in Ref. p2f . 

In principle, the LEC's of Eq. ( |C1|) determine, or viceversa they can be determined from, 



the next-to- leading order results of ChPT, as it is explicitly shown in Ref. for the case 
of elastic tttt— scattering. However, the perturbative calculation is not available and it only 
exists next-to-leading results for ttN, with no coupled channels |],|I7],[^,[?(J , in the framework 
of Heavy Baryon Chiral Perturbation Theory (HBChPT). In Appendix [D] we discuss this 
point in more detail. 



C. Predictions for other processes 

In Fig. |] we show some of our predictions for phase shifts, inelasticities and s— wave 
T = 1/2 partial cross sections for some other channels. For most of them there are no 
data. The 7]N — > rjN elastic phase shifts (top left pannel) present a steep raise close to the 
threshold going up to 70°, which corresponds to a typical low energy resonance behaviour 
triggered by the iV(1535) resonance (see next section) . Accordingly, the corresponding 
partial cross section (bottom pannel) takes a unnaturally large value as compared to other 
elastic and transition cross sections. This is in contrast to any expectation based in the Born 
approximation, since the corresponding potential in this channel vanishes (see Eq.( ^)). 
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FIG. 4. Top pannel: s— wave T — 1/2 phase shifts (in degrees) for elastic ijN — > rjN (left pannel) and 
KA — » KA (right pannel) processes as functions of the CM energy. Middle pannel: same as before but for 
inelasticities. Vertical lines indicate the opening of reaction thresholds. Bottom pannel: s— wave T — 1/2 
meson-baryon cross sections (7r7V — ► 7r7V, 77 iV, KA, KT,; rjN — > ryiV, KA, KE; KA — » KA, KS) in mbarns as 
functions of the CM energy. Dashed lines indicate nN initial state. Solid lines indicate rjN initial state. 
Crosses indicate KA initial state. All the lines start at the relevant final state threshold (with the exception 
of the elastic irN — > nN reaction). 
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The effect of the iV(1535) can also be seen at the figure (bottom pannel) by the maximum 
in ttN — > r]N cross section and the cusp effect in the ttN — > ttN partial cross section. On the 
other hand, KA — > KA phase shifts (top right pannel) turn out to be extremely small. The 
effect of the iV(1650) can also be seen at the cross sections, particularly in the 7riV — > ttN 
partial cross section, though the effect is less pronounced than in the iV(1535) case. We 
have not plotted the elastic ICE —>■ KT, cross section since the physical process involves also 
the isospin T = 3/2 channel, not considered in this work. 

Our estimates for r)N and KA scattering lengths (defined for the elastic channels similarly 
as in Eq. ( [5gD are 

a vN = 0.772(5) + i 0.217(3) fm 

a KA = 0.0547(5) + i 0.032(4) fm (39) 

respectively. The scattering length a V N compares reasonably well with the one obtained in 
Ref. @, <Vv = 0.68 + i0.24 fm. 



D. Second Riemann sheet: poles and resonances. 

In this section we are interested in describing masses and widths of the Sn— resonances. 
An illustrative picture of the complex CM energy plane with the singularities from the 
Particle Data Book pTI is presented in Fig. |5]. For a more distinctive characterization of the 
resonances one has to look for poles in the complex s— plane. 



Since causality imposes the absence of poles in the physical sheet |[28|| , one should search 
for complex poles in unphysical ones. Among all of them, those closest to the physical 
sheet are the most relevant ones. For the sake of clarity, we will devote some space here 
to explain, in a quantitative manner, the meaning of "close" in this context. We look for 
poles in the coupled channel matrix amplitude t(s) defined in Eq. (|3"3"D. We have only 
examined the entry 11, ttN — > ttN, of that matrix. The position of the complex poles, as 
long as they are produced for physical resonances, should be independent of the particular 
channel. However, residues at the pole do depend on the examined channel, because they 
determine the coupling of each of the channels to the resonances. This interesting point will 



be discussed elsewhere 29 



To begin with, let us assume a situation where the coupled channel formalism is not 
needed, i.e. an artificial situation where only the 11 element of the first column and row of 
the D— matrix in Eq. (|8]) is non-vanishing. In such a case, elastic unitarity requires only a 
unique finite branch point at s = s + = (M N + m n ) 2 and a cut along the line [s + , +oo[. The 
scattering amplitude in the unphysical second Riemann sheet (i/j(s)) is simply obtained by 
analytical continuation of the amplitude in the physical first Riemann sheet (t(s) = i/(s)) 
across the unitarity cut, and therefore the following relation for inverse amplitudes should 
hold (s real and above s + ) 

tj 1 1 {s + ie) = tj 1 {s-ie) (40) 



18 



Complex CM Energy Plane 



7tN 



T|N 



KA 



KZ 



x 

N(1535) 



N(1650) 



FIG. 5. Location of reaction thresholds and resonances in the complex CM energy plane. The corre- 
sponding unitarity cuts have increasing thickness for increasing energy 



The unitarity condition for the inverse amplitude, deduced from Eqs. 
reads, 



Disc [t _1 (s)] = *7 1 (s + ie) - tj 1 ^ - ie) = 2ip(s) 
with s G R, where the phase space function 



s > s + 



X l / 2 {s,m w ,M N ) {^~s + M N f -ml 

— X 



167TS 



2v^ 



and (PflSD 

(41) 

(42) 



has been introduced, understanding that p(s) is a function of the real variable s. Then, 
analytically continuing the phase space function to all complex plane, the unitarity condition 
reads 



t7 x (s + ie) - tj 1 ^ - ie) = 2ip(s + ie) s + < s G 11 



(43) 



where the cuts for \?(z, m%, M%) go along the real axis for — oo < s < s_ and s + < s < oo. 
The function is chosen to be real and positive on the upper lip of the second cut, s + < s < oo 
and hence it satisfies: 



\'(s + ie,n4,M N/ 



\\2(s,ml,M^j 



s + < s G R 



(44) 



Besides, the cut for the ^fz function, also appearing in p(z), is taken along the line ] — oo, 0] 
and the multivalued function is taken to be positive for real and positive values of z. Now 
using Eqs. (flOP and ( |43|) one finds the amplitude in the second Riemann sheet, 



tjl(z) = tj\z) - 2ip(z) z G C 



(45) 



On the other hand, 
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tjj(s - ie) = tj 1 ^ - ie) - 2ip(s - ie) =t I 1 (s + ie) - 2ip(s + ie) - 2ip(s - ie) 

= tj 1 (s + ie) s + < s E R (46) 

which means that there are only two Riemann sheets linked to the unitarity cut. The 
analytical structure, concerning the unitarity cut, of the inverse amplitude is determined by 
the function 

- Ms) x C~ m ^ s Ml + M N ) (47) 



as it is deduced from Eqs. fl33|), fl30|) and ( |A8| ) . Because of the choice of cuts for the 
multivalued function \fz above, the function L(s) introduced in Eq. ( |A9|) determines the 
analytical structure of t(s). Thus, the two Riemann sheets of t(s) related to the unitarity 
cut, are obtained from the values n = and 1 in Eq. (|A13| ) for L(z,n). 



In the general case of multiple thresholds, as it is the case in this work, the above 
conclusions hold for any of the finite branch points, located at the physical thresholds, and 
unitarity related cuts going along the lines [branch point, oo[. 

For four channels, there are a total number of 2 4 = 16 Riemann sheets related to unitarity. 
Now, Jq(s) and hence L(s) become diagonal matrices in the coupled channel space, namely 



L(z,n) 



/ ' L m „M N (z,n nN ) \ 

£m 1) M J v(^,?Vv) 

L mKMA (z,n KA ) 

V L mK M s {z,n K s) J 



(48) 



with z G (D, n = (n n N, n V N, uka, tlkt.) and the dependence in the masses of the function L 
is explicitly given. The first Riemann sheet (i/(s)) corresponds to the choice n = (0, 0, 0, 0). 
As mentioned above, poles can only occur in any of the remaining 15 Riemann sheets. The 
closer the position of the pole to the scattering region (upper lip of the first Riemann sheet) 
the bigger is the influence on the scattering amplitude. All Riemann sheets can be reached 
continuously from the first one by looping around the appropriate branch points. From 
this point of view, close means proximity following a continuous path. Thus, for the region 
(m n + Mm) 2 < s < (m v + M^) 2 the poles of the (1, 0, 0, 0)— Riemann sheet located in the 
fourth quadrant and with values of Rez belonging to the above interval, are expected to have 
the biggest influence on the scattering amplitude, as it is illustrated in Fig. |6|. Similarly, for 
the region (m v + M N ) 2 < s < {m K + Ma) 2 the poles of (1, 1, 0, 0)— Riemann sheet located 
in the fourth quadrant and with values of Rez belonging to the above interval, are expected 
to play a crucial role, and so on... 
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Im(z) 




FIG. 6. Different paths, in the s— complex plane, showing how to reach the point S, located in the 
physical scattering region in the interval [(m ff + Mm) 2 , {m n + Mjv) 2 ] , from points P (eventually poles) 
located in different Riemann sheets, denoted by the vector n as introduced in Eq. (J48|), and placed both in 
the first and fourth quadrants. The unitarity cuts are also depicted in the figure. The "distance" between 
S and P is obtained by the length of the shortest path joining them. This can be achieved after continuous 
deformation of the paths depicted in the figure, i.e. any deformations without intersecting the branch points. 

Thus, we define the " Second Riemann Sheet" in the relevant fourth quadrant (tri(s)) 
as that which is obtained by continuity across each of the four unitarity cuts. It is obtained 
using for the diagonal matrix Jq(s) the following function^ 




13 Though, each of the functions L(z;n) are analytical in the complex plane, except for the per- 
tinent unitarity cuts, note that Cjj so defined, is continuous for real values of s, but presents 
additional discontinuities out of the real axis. 
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FIG. 7. Modulus of the ttN — ► ttN element of the scattering amplitude t(s), defined in Eq. (|33|), in 
the s— complex plane. In both plots, vertical lines indicate the position of the poles. Top panel: Fourth 
quadrant of the "Second Riemann Sheet", as defined in Eq. (ff9|), and the physical scattering line. The two 
observed poles are identified to be the Six— iV(1535) and — iV(1650) resonances as it is discussed in the main 
text. Bottom panel: Fourth quadrant of the "Second Riemann Sheet" and the first quadrant of the first 
(physical) Riemann sheet. Besides the two poles already appearing in the top panel, there is a third one. 
Though it is unphysical because it appears in the physical sheet out of the real axis, it does not influence 
the scattering line as the plot clearly shows. 



22 



In Fig. [7] we show the absolute value of the tn(s) element of the scattering amplitude 
both for the fourth quadrant of the "Second Riemann Sheet" and the first quadrant of the 
first (physical) Riemann sheet. The physical scattering takes place in the scattering line in 
the plots (upper lip of unitarity cut of the first Riemann sheet). The positions of the two 
poles in the "Second Riemann Sheet" are (s = M R — iMrTr): 

First Pole : M R = 1496.5 ± 0.4 T R = 83.3 ± 0.7 (50) 



Second Pole : M R = 1684.3 ± 0.7 T R = 194.3 ± 0.8 (51) 

where all units are given in MeV and errors have been transported from those in the best 
fit parameters (Eq. (pT|)), taking into account the existing statistical correlations, through 
a Monte-Carlo simulation. These poles are resonances and can be identified to be the Sn — 
iV(1535) and — iV(1650) ones which, according to Ref. [^] (PDG), are located at 

iV(1535) : M R = 1505 ±10 T R = 170 ± 80 (52) 



iV(1650) : M R = 1660 ±20 T R = 160 ± 10 (53) 

where again units are in MeV and we quote data from position of the poles which are 
slightly different to those deduced from a Breit-Wigner fit. The agreement of our predic- 
tions and the PDG ones is satisfactory. Our calculated width for the iV(1535) / iV(1650) 
turns out to be smaller /larger than the experimental one, in great part, because its mass 
is slightly smaller /larger than the data and hence the available phase space for the decay 
decreases/increases. Besides, the inclusion of the three body channel (tttcN) would influence 
both the widths and the masses of the resonances. 

Residues at the poles depend on the examined channel, because they determine the 
coupling of each of the channels to the resonances. Thus, from the results shown in Fig. (^) 
one could predict the coupling of the iV(1535) and iV(1650) resonances to nN. A detailed 
study of the couplings of these resonances to all channels, not only nN, is presently under 
way. 

On the other hand, there is a unphysical pole in the physical (first) Riemann sheet. 
It is located at (s = M 2 + IMF) M « 1582 MeV and T « 166 MeV and it violates the 
Mandelstam's hypothesis of maximal analyticity P5 |. This unphysical pole appears because 



we have truncated the iterated potential to solve the BSE. However, as can be seen in the 
plots, the two poles in "Second Riemann Sheet" have a much larger influence on the physical 
scattering than the spurious (unphysical) one. Thus, the influence of this unphysical pole 
may be disregarded^. 



Due to Schwartz's Reflection Principle there is also a pole in the fourth quadrant of the first 
Riemann sheet (complex conjugated of that given above) which influence is even more negligible 
than that of the first quadrant. This is because it is placed at a substantially larger distance of 
the upper lip of the unitarity cut. Existence of other complex conjugated poles, both in the first 
Riemann sheet or in any of the unphysical sheets, is not precluded, but from Fig. (R) we infer that 
their influence in the scattering is not significant. 
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IV. CONCLUSIONS 



In this paper we have developed a Bethe-Salpeter formalism to study s— wave and T = 
1/2 meson-baryon scattering up to almost 2 GeV. We work on a four dimensional two body 
channel space and the kernel of the BSE takes into account chiral symmetry constraints as 
deduced from the corresponding effective Lagrangian. At lowest order in the chiral expan- 
sion for the potential, an analytical explicit solution is found which manifestly complies with 
multiple channel unitarity. Among the several issues which can be explored using the present 
formalism, we have focused our attention on nN elastic scattering (phase shifts and inelas- 
ticities), and the measured inelastic cross sections, being the agreement with experiment 
rather good. Besides, some predictions for other cross sections, not yet measured, have been 
also given. We have undertaken a careful discussion on the analytical structure and continu- 
ation of the scattering matrix amplitude to the complex s— plane, which becomes mandatory 
in order to extract the location of the Sn resonances. We have searched for poles in the 
"Second Riemann Sheet" and compared both masses and widths to data. The agreement 
is also quite satisfactory. Thus, and despite of having neglected the three body production 
channel tcttN, we provide a rather successful description of the s— wave and T = 1/2 meson 
(it, rj, i^)-baryon (N,A, S) scattering up to almost 2 GeV in the strangeness zero chan- 
nel. Couplings of the iV(1535) and iV(1650) resonances to the different open meson-baryon 
channels can be obtained from our amplitudes and it will be discussed elsewhere. 

Nevertheless our calculation has some obvious limitations and hence it might be im- 
proved. Besides the inclusion of the itnN channel^, one should consider the inclusion of 
higher order terms in the two particle irreducible matrix amplitude (potential) which would 
lead to a realistic predictions for higher partial waves in the T = 1/2 channel. Understand- 
ing the free parameters of our model and their numerical values presented in Eq. ( |Cl~l) , in 
terms of the LEC's appearing in the higher chiral Lagrangian order pieces, either within the 
HBChPT formalism [17] or in the fully covariant framework recently developed in Ref. [|30 



would be obviously desirable. As we have shown in Appendix |B|, such a task would, at 
least, require to know the perturbative l// 6 order to take into account the likely important 
rjN— channel effects, which are effectively incorporated in the BSE scheme presented in this 
work. 

Given the phenomenological success of the presented framework, it seems natural to 
extend it to other non zero strangeness channels, for a recent overview of related aspects see 
for instance Ref. fl3l| , or to the study of meson photo-production processes. 



In our point of view, it is a highly non trivial task to find a solution of the BSE including 
a three body intermediate state exactly complying to three body unitarity. In some cases, for 
instance for elastic ttN — > ttN, meson rj production, some insight might be obtained by treating 
perturbatively the process, as our results suggest. 
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APPENDIX A: BASIC INTEGRALS 

We display the explicit expressions for the loop integrals used in this paper. The basic 
integrals appearing in the solution of the BSE are: 



J ^^!wr¥^-J^M (A1) 
J *<« = '/ (§v^ jzjzz* = - *> - a,,, (A2) 

J ^ = '/ (Sv^^j^ = {r ~ *' wmr ~ *> ~ {r ~ * )A * (A4) 

and the results are obtained from relativistic and translational invariance requirement in 
momentum space. Here, is a quadratically divergent integral 

= (S)V^ <A5) 

which would require renormalization. Besides, the linearly divergent integral J{f) can be 
evaluated yielding: 



j(f) = r 



2s K ' 2s 



+ MJo(s) (A6) 



A A A - ' f ^ q 1 . f d 4 q 1 

where is quadratically divergent as well and the logarithmically divergent integral 

Jo(s) needs one subtraction to make it finite. Choosing for definiteness the threshold value 

s = (rh + M) 2 we get 

with J (s — (m + M) 2 ) a divergent integral and the finite function Jq(s) is given by0 



16 Jq(s) is a diagonal matrix in the coupled channel space and for simplicity we work from now on 
in a given channel. 
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Ms) 



(4tt) 



M 2 — m 2 M -m 



M + m 



In— + L(s) 
m 



(A9) 



and for real s and above threshold, (m + M) 2 , we have 



L(s) = Lis + ie) 



\ l/2 (s,m 2 ,M 2 ) 



log 



S — S- 



(A10) 



where have defined the pseudothreshold and threshold variables as 

s_ = (m - M) 2 s + = (m + M) 2 



(All) 



respectively, and the logarithm is taken to be real. Note that L(s + ) = 0. For s > s + the 
imaginary part along the unitarity cut may be computed directly from the above Eq. (|A10|) 
or through Cutkosky's rules, 



2i ImJ (s 



Disc Jo (s) 
d A q 



-2i 



(2tt) 2 

. A 1 / 2 (s,m 2 ,M 2 ) 



[J (s + ie) - J (s - ie)] 
2m) 2 5 + (q 2 - m 2 )5 + ((P - q) 2 - m 2 ) 



167TS 



(A12) 



Up to a 167r 2 factor the function L(s) has the same discontinuity as the function Jo(s). 
Taking into account that we have to evaluate the function Jq(s) not only for real s > s + but 
also below threshold^] and in the second Riemann sheet as well, to look for the position of 
resonances in the complex s— plane, we give here the analytical continuation of L(z) used in 
our calculation. Defining p± = \z — s±\ and taking the principal arguments, Arg(...), 9± of 
\z — s±\ to lie in the range < 9 + < 2n and — 7r < 6- < n respectively we have 



L(z, n) 



R{z) 



{p+P- 



>l/2 



,i(6» + +6»_+2n 



*V 2 {hi\R{z)\+iAxg[R{z)\-2m) 



p V2 e W + /2 + p l/2 e i0_/2 e m 7 r 
p l / 2 e i6 + /2 _ p 1 / 2 e i6»_/2 e in7r 



(A13) 



where Arg[R(z)} should be taken in the interval [0, 2tc[. For n = one gets the first Riemann 
sheet Lj(z) = L(z,n = 0), which only has a (unitarity) cut along the real axis s + < s < 
oo. When going across the unitarity cut once we jump into the second Riemann sheet, 
corresponding to n = 1, Lu(z) = L(z,n = 1). If we loop twice around the threshold branch 
point z = s + we come back to the original Riemann sheet. The second Riemann sheet has 
an additional cut along the real axis — oo < s < s_ and the following relation holds 



17 For instance, when calculating the elastic irN scattering, obviously there are values of s below 
heavier thresholds r]N, KA, KT,. 
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L n {z) = L^z) + 2?ri — ^ 2 i (A14) 

z 

where the cuts for A3 (z, m 2 , M 2 ) go along the real axis for — oo < s < s_ and s + < s < oo. 
The function is chosen to be real and positive on the upper lip of the second cut, s + < s < oo 
and corresponds \\s (s, m 2 , M 2 )\ = A2 (s + ie, m 2 , M 2 ). The — 2m constant appearing in 
Eq. ( |A13| ) determines the chosen Riemann sheet of the logarithm and ensures that Lj(z) 
is purely real along the real axis below threshold. Note that since R(z) only vanishes at 
infinity we never have a chance to cross the cut of the logarithm and never change log- 
Riemann sheets. 



APPENDIX B: DERIVATION OF THE SOLUTION OF THE BSE 

Here we show how to derive Eq. (f28[) displayed in the main text. The ansatz of Eq. ( |26|) 
reduces the BSE integral equation Eq. ( |20| ) into a set of linear equations for the matrix 
coefficients a,b R ,bL and c 

a = a[(f - M)J - A^ + b L (f - M)[(f - M)J - A^]-^ 
b R = b R [{f - M)J - A^ + c(f - M)[(#> - M)J - A A ]^ + j 2 
b L = aJj 2 + b L [{f - M)J - A^} — + j 2 



c = b R Jj 2 + c[(f - M)J - A„]^ (Bl) 



In the above equation J stands for J(f ) defined and evaluated in Appendix 0. The solution 
of this matrix system is tricky although straightforward. The main complication arises 
from the non-commuting character of the fermion mass matrix M with the coupled channel 
matrices a, &l,6r, c and D. Defining 



X = a Y R =(f-M)b R Y L = b L (f-M) 

1 

1 

The set of matrix equations can be written as 



Z=(f -M)c(f-M) G = — (f-M)J-A^ (B2) 



X = XGD + Y L GD (B3) 
Y R = Y R GD + ZGD + ^ ~ f )D (B4) 

r 

Y L = X(G + ^)—^D{f — M) + Y L G—^-D(f — M) + °^~ M) (B5) 
P f-M f-M P 

Z = Y R (G + ^) — l -^-D(f - M) + ZG — l -^D{f-M) (B6) 
P f-M f-M 



Summing Eqs. ( p3|) and (|B4] ), and Eqs. (|B5| ) and (|B6| ) , we get after some matrix reshuffling 
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(f — M) 
P 

(f-M) 



(X + Y R )(D- 1 -G)-(Y L + Z)G 



r , (X + Y R )(G + ^)-(Y L + Z) 



M)~ 1 D~ 1 ' 



M)-G 



(B7) 
(B8) 



Subtracting and summing Eqs. 



and 



wc 



(X + Y R ) = (Y L + Z) 



1 n-l, 



M)~ D 



(X + Y R )(D- 1 -2G 



A T 



M)(p- X 



+ (Y L + Z) 



1 m \ — l 



P 



l n-l/ 



M)~ D 



P y nJy ' P 

respectively. We can then solve for Z + Yl from Eqs. ( |B9|) and (|B10|) yielding 

A 



(B9) 

M) - 2G] (BIO) 



(n + ^r 1 = { 



Mp L D- L (f - M)(D 



-1 i ^niN-i 



./ 2 



-2G) 



+ (f — M) D (f — M) 



P" 1 



A, 



./ 2 



A r 



/ 2 



1 



/ 5 



2(f - M) 



(Bll) 



Using the proportionality relation between X + Y# and Yl + Z given by Eq. ( |B9| ) we 
obtain the following expression for the on-shell t(f ) matrix, 



t{f ) = x + Fr + y L + z = (y l + z) 



A 



m \ — l 



/ 2 



1 



(B12) 



Inverting this equation and using Eq. ( p 1 1|) we get finally the expression given in Eq. (f2~8|). 



APPENDIX C: BEST FIT RESULTS 



The best fit (x 2 /dof = 0.75) parameters are 



J; 



A A r i i/(m 7r 

A A /(m x 
A s /(m K 
A 7r /(m 7r 
A ri /(m v 
A K)l /(m K 



7r7V 

Jr]N 

Jka 
Jkt, 



+ M N f 
+ M N f 
+ M A ) 2 
+ M S ) 2 
+ M N f 
+ M N f 
+ M A ) 2 
+ M 2 ) 2 



0.1897 ±0.0004 
0.6206 ± 0.0002 
-1.227 ±0.003 
-0.0143 ±0.005 
0.776 ± 0.002 
1.8375 ±0.0004 
-2.923 ±0.008 
1.000 ±0.017 
-0.0123 ±0.0003 
-0.1560 ±0.0002 
-0.006324 ± 0.000003 
0.001128 ±0.000003 



(CI) 
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The correlation matrix, defined as usual, 



.1 ; 



(/(Ci, . . . , C n )) 



(Ci) 



\[&) - <^> 2 

1 N 

T7 /( C 1,C5 • • • > c n,a) 

^ Q=l 



(C2) 



being q any of the 12 parameters: J's and A's, turns out to be: 
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1.00 
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1.00 
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0.21 


1.00 
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-0.28 


0.48 


1.00 
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-0.22 


0.57 


0.09 


1.00 
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-0.35 


-0.18 


-0.15 


0.00 


1.00 












0.32 


0.24 


0.97 


0.42 


0.49 


-0.10 


1.00 










0.36 


-0.30 


0.51 


0.99 


0.14 


-0.09 


0.43 


1.00 
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-0.30 


-0.08 


-0.21 


0.91 


-0.25 


-0.04 


1.00 






-0.51 


0.48 


-0.48 


-0.62 


-0.71 


0.30 


-0.33 


-0.65 


0.30 


1.00 
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0.10 


-0.37 


0.12 


0.20 


0.56 


-0.59 


0.01 


0.18 


-0.55 


-0.77 


0.19 


-0.03 


0.54 


0.48 


0.67 


-0.56 


0.42 


0.48 


-0.69 


-0.84 



1.00 
0.72 



(C3) 



1.00 



The large correlations (0.97 and 0.99) between the pairs Jka — A A and Jky, — A^ can be 
understood by looking at Eq. ([34]) and taking into account the smallness of the parameters 
Ajfi and Ak,2 respectively. 



APPENDIX D: THE HEAVY BARYON LIMIT AND HBCHPT 

It is in principle of theoretical and phenomenological interest the study of the heavy 
baryon limit of the BSE amplitude given in Eq. fl34|). In the static limit, baryons behave 
like fixed sources, and consequently the two particle problem should reduce to a one particle 
scattering problem (in our case of meson-baryon scattering it would correspond to a Klein- 
Gordon equation with a spin-dependent potential). It has been known that the BSE has some 
difficulties in reproducing this heavy-light limit in certain situations (ladder approximation 
to one boson exchange |32[]). We show below that our amplitude has a correct one particle 
limit, due to the fact the s— wave interaction is of the contact type. If, in addition to a heavy 
baryon expansion, a chiral expansion in powers of 1/f 2 is carried out, we should recover in 



this double expansion some form of the results found in Ref. p5| , |26| , ^4| within HBChPT for 
the elastic irN scattering amplitude. The heavy baryon expansion may be taken by making 
the baryon masses M — > oo but keeping the meson masses, m, and the meson momentum, 
q, finite. On the other hand, baryon mass splittings must be considered higher order effects, 
see e.g. [19|1, so that we take the mass matrix 

M = M B + AM (Dl) 

with Mb — > oo the common mass of the baryon octect which is proportional to the identity 
matrix. Accordingly, in the ttN elastic channel we take 
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yfs = E + uj = M n + uj+ * + • ■ ■ (D2) 

where Mn = Mb + AMn- In this appendix we match our amplitude to the HBChPT third 
order results of Ref . [TJJ] based on the previous analysis of Refs. [^,0] . The heavy baryon 
expansion can be directly done for explicit variables, such as masses and CM energy ^fs. 
The constants J-^n, • • • , Ajv,i, • • • , A w , . . . , given by our numerical fit in Eq. ( |Cl|) , also might 
have a baryon mass dependence. Such a dependence should lead to some changes in the 
heavy baryon expansion which cannot be easily determined. In addition, given the non- 
perturbative nature of our solution and the fact that many aspects of the coupled channel 
meson-baryon data in the Sn-strangeness zero channel are described after unitarization up 
to energies as high as -/s = 2GeV, it seems obvious that the parameters of Eq. (|C1|) also 
incorporate higher order effects in the chiral expansion. 



1. Static limit 

It is convenient to do the study in terms of the inverse coupled channel matrix ampli- 
tude, t _1 (s), given by Eq. fl3~4|). From the expression of the one-loop integral Eq. (|A9|) and 
Eq. (|A10|) we get to leading order 

1 ( M 2 \ - / 1 \ 

MJ (s, m, M)l^ V Mj ^ 2 _ m ^ = — log j (m - u>) + K m {u) + O [-) (D3) 

The subtraction point for the HBChPT integrals is uj = m, K m (m) = 0, and their explicit 
expression is 



'uj 2 


— m 2 


bo 


-UJ 2 



-m 



K m (uj) = — I +Vm 2 -w 2 arccos (-^) uj 2 < m 2 (D4) 



uj > m 

Thus, one obtains for the unsubtracted integral 

MJ ( g ,m,M)|^ =V M ^ 2 _ m ^ = K mM (uj) + o{^ 

1 M 2 - / I \ 

= MJ ™m + T^-2 h S^(m-uj) + K m {uj) + O [-) (D5) 

where the HBChPT unsubtracted integrals K m M(uj), fulfilling K m Mijn) = MJ m u and the 
heavy baryon approximation of the subtraction constant defined through Eq. 



J (s = (m + M) 2 ) = J° mM {l + O } (D6) 



have been introduced. In the static limit we obtain from Eqs. (16) and (|34]) (f(oj) 

-t(s)/M ) 
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47T 



167T 2 



M 2 
In — — (m 



u) + M.r, + 



A - 
4M 



A u ~ - 

m 



P 



:D 



f 



DA°~D 

4 m 



-1 



(D7) 



with Kfaiuj) = Diag(K n (uj) , K v (u), Kk{oj), Kk(uj)) and the heavy baryon approximation of 
the subtraction constants are defined by means of the expansion 

A, fl = <{l + (±)} 

The Eq. ( p7|) corresponds, as it should, to a one particle scattering problem, fulfilling the 
coupled channel unitarity condition 



Im/(co>) 1 = —\Juj 2 — m 2 Q(uj — 



m) 



(D8) 



The pole in Eq. (|D7| ) for the inverse amplitude is a static limit reminiscent from the baryonic 
Adler zero, y/s — M = 0, of the lowest order potential. The constant combination appearing 
in the inverse amplitude, Eq. ( P7|) , M J°~ fa + A°^./4M should go to some definite value in 
the static limit, M — ► oo. In case it would diverge, the scattering amplitude would become 
trivial. We may try to estimate these constants using the numerical values obtained in the 
X 2 — fit carried out in this work and given in Eq. (|CT|). We get 



M A J KA + 



A^tv 


= -0.47m 


AM N 




A^jv 


= -1.08m 


AM N 




A^a 


= +0.67m 


4M A 


A^ s 


= -1.24m 


4M S 





(D9) 



in units of the relevant pseudoscalar meson masses. Though the numerical values used for 
the subtraction constants contain higher order effects in the heavy baryon expansion, we see 
that there is indeed some trend to cancellation, because the J's and the A's contributions 
have opposite signs in the first three cases, and Jkt, is very small. Moreover, the constants 
do not seem to attain unnaturally large values, although it is hard to say which should be 
an accurate appropriate value. 



2. Chiral and heavy baryon expansion 

Expanding Eq. fl3~4T) in powers of 1/ f 2 we get 

t(s) = t 2 (s) + Us) + ■ ■ ■ 



(D10) 
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where 



1 r 
p{ 


\/i — M, -D} 


-if 


v/i M) D- 


1 r 

71! 


Vi-M,£>} 



M 



\sfs + M) 2 - m 2 

, V5 



Jo(s) 



A, 



+ 



*rhM 



(Dll) 
(D12) 

M,d} 



In the heavy baryon limit we get for the elastic nN - 

f(co) = f 2 {u) + U(u) + 



y/s-M ' 2Vs , 
ttN amplitude in the Su channel 

(D13) 
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(D14) 
A^ " 



To~Ak-, 2 
/ 



+ 



2M E , j 
(D15) 



In the regionQ < a; < m# only the K n N(co) has an imaginary part, to comply with 
perturbative elastic unitarity, whereas Kka_{oj) and Kkh{oj) are purely real. To write down 
this expression we have considered the prescription D/f 2 ^ f~ l Df~ x given in Eq. (|9]). At 
threshold, uo = m n , the ttN scattering length in this channel reads 
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(D16) 



From Eq. ([13 1 5| ) and Eq. ( |D16|) above it is clear that there is no contribution from the rjN 
channel to this l// 4 order of approximation. This is a direct consequence of the structure of 
the coupled channel matrix D given by Eq. (§) since the corresponding ttN — > rjN transition 
matrix element vanishes in the Born approximation. For this reason, the channel rjN starts 
contributing at order l// 6 . This situation is unexpected, because on general grounds we 
expect the r]N channel to be more important at low energies than the KA and KH channels 



18 Note, that the r/N channel appears at order l// 6 . 
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since their thresholds lie at higher energies. Using the numerical values of the coefficients 
obtained from the fit of Eq. flClQ to estimate the scattering length, Eq. (P16 ), we get 



nN KA KS 

a wN = 0.22 + 022- L06 + 018 = -0.43 fm (D17) 

V/ 2 i// 4 

which should be compared to the one of Eq. (|3|) obtained from the full amplitude, 0.179 fm. 
Obviously, the different values should be attributed to non negligible higher order effects, 
which in particular include rjN contributions, and higher order corrections to the J° and 
A coefficients. This is very reassuring because the full BSE amplitude, Eq. (|34"D, besides 
restoring unitarity automatically includes all orders in the chiral expansion. 



3. Matching to HBChPT 



The 
to third 

/a(w) 



Sii partial wave amplitude deduced from the work of Ref. based in HBChPT 
order [25|,26|j reads, after straightforward angular integration, 
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(D18) 



(D19) 



There is no unique way to match the low energy chiral expansion of the coupled channel BSE 
amplitude, Eq. ( [Djg) , to the third order HBChPT calculation of Refs. ^,H, Eq. ( pT9|) . 
The analytical structure is different besides the elastic unitarity cut at uj = m n which 
turns out to coincide. Indeed, while the former presents the inelastic unitarity cuts for the 
considered KA and KT, coupled channels, the latter includes perturbatively the left hand 
cut at uj = 0. Obviously, any particular choice of the matching point generates a specific set 
of low energy constants. After explicitly separating the elastic unitarity correction of both 
amplitudes, it seems reasonable to do the matching of the remaining pieces in a polynomial 



expansion around the elastic threshold point, tu 



since neither inelastic unitarity 
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cuts nor the left cut are expected to be crucial at that point. Instead, we expect both 
amplitudes to provide a sensible approximation. Also, direct inspection of Eq. ( P15| ) and 
Eq. ( pi9| ) reveals that only some additive combinations among renormalization constants 
can be established. In particular, in Eq. (|D15|) there are two independent combinations of 
low energy constants. Thus, it proves sufficient to Taylor expand around u = m n up to first 
order. Using the numerical values for the input parameters, Eq. (|37|) , in Eq. (|D15|) and the 
numerical values for the parameters in Ref. [|nj] in Eq. ( pi9| ) the following identifications 
hold, (in units of fm) 

0.185 = 0.175 + 22.8A° + 2.51A^ 1 + 17.3A^ 1 + 1.57A A + 2.12A^ 2 + 0.179A£ 

- 7MJ° N - 3MJ° KA - 0.357J£ S 

-0.051 = 0.163 + 20.3A° + 5MA° N1 + 15.7 A° K1 + 3.14A A + 1.94A^ 2 + 0.359A| 

- 15.2J° N - 6.01J^ A - 0.714j£ s (D20) 

where A^ = A%/(m + Mf and = A^/(m + M) 2 are dimensionless. As we see, there 
is a large degree of redundancy when the matching is performed considering only these low 
orders of the expansion. By using the values of Eq. ( |C1|) , to estimate the heavy-baryon mass 
independent parameters appearing in the right hand side of the two relations established in 
Eq. (|D20|) , we obtain —0.43 and —0.618 respectively. The disagreement, with respect to the 
left hand side values, is not completely surprising because the numerical values used for the 
subtraction constants contain higher order effects in the heavy baryon expansion. Besides, 
the nominally small differences (0(1/M)) between A°'s and J°'s and A's and J's might lead 
to significant numerical changes because the factors multiplying these constants are large in 
units of the left hand side values of Eq. ( |D20| ) . 
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